Diffusive Growth of a Single Droplet with 
Three Different Boundary Conditions 

Z. Tavassoli and G. J. Rodgers 

Department of Mathematical Sciences, Brunei University, 
Uxbridge, Middlesex, UB8 3PH, U.K. 

Abstract 

We study a single, motionless three-dimensional droplet growing by ad- 
sorption of diffusing monomers on a 2D substrate. The diffusing monomers 
are adsorbed at the aggregate perimeter of the droplet with different bound- 
ary conditions. Models with both an adsorption boundary condition and 
a radiation boundary condition, as well as a phenomenological model, are 
considered and solved in a quasistatic approximation. The latter two models 
allow particle detachment. In the short time limit, the droplet radius grows 
as a power of the time with exponents of 1/4, 1/2 and 3/4 for the mod- 
els with adsorption, radiation and phenomenological boundary conditions, 
respectively. In the long time limit a universal growth rate as [t/ ln(t)] 1//3 
is observed for the radius of the droplet for all models independent of the 
boundary conditions. This asymptotic behaviour was obtained by Krapivsky 
||T5fl where a similarity variable approach was used to treat the growth of a 
droplet with an adsorption boundary condition based on a quasistatic ap- 
proximation. Another boundary condition with a constant flux of monomers 
at the aggregate perimeter is also examined. The results exhibit a power law 
growth rate with an exponent of 1/3 for all times. 
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1 Introduction 



There has been considerable interest in diffusive-growth processes including 
growth phenomena for a droplet on a substrate. This growth phenomena in 
the case where diffusion and coalescence play the major roles are common 
in many areas of science and technology [|IJ, 0. In this process, each droplet 
diffuses and grows individually and coalesces with contacting droplets. The 
kinetics of these phenomenon have been studied experimentally and theoreti- 
cally [3-29]. Some models have been developed to explain the kinetics of these 
processes. One such model [§, [12], [y| consists of a single, motionless three 
dimensional droplet formed by diffusion and adsorption of non-coalescing 
monomers on a 2D substrate. In the model it is assumed that the diffusing 
monomers coalesce only with large immobile growing droplet and not with 
each other. In || [12| a static approximation was used to solve the diffusion 
equation and an approximate description of the long time behaviour was 
obtained. The static approach predicted an asymptotic power law growth 
rate for the radius of the droplet. Because of the growth of the droplet, the 
present problem involves a moving boundary. Moving boundary problems 
in the context of the diffusion equation are referred to as Stefan problems 
[30-32]. The only exact solutions for these problems have been found using 
a similarity variable method, see for instance, [30-34] and references therein. 
Using this method for a droplet of dimensionality d growing on a substrate 



of the same dimensionality, an exact scaling solution can be found. In |33 



such a solution in one dimension has been derived which can be generalised 
to a higher dimension. However, the problem of a 3D droplet growing on a 
2D substrate, may be treated by approximate methods. A simple treatment 
based on a quasistatic approximation has been presented in ||15||. A simi- 
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larity variable approach was used to solve the Stefan problem with moving 
boundary |TJj . The results predicted that the radius of the droplet increases 
as [t/ baft)] 1 / 3 asymptotically. The asymptotic growth law predicted by the 
static approach differs from the quasistatic answer by a slowly varying log- 
arithmic factor. In all the models in [El [T2|, |TjJ an adsorption boundary 
condition at the aggregate perimeter of the droplet, was considered. 



In [37 a generalisation of Smoluchowski model |36| for diffusional 
growth of colloids, was presented. Smoluchowski |35| considered the process 
of diffusional capture of particles assuming the growing aggregate is mod- 
eled as a sphere. He then solved the diffusion equation with an absorbing 
boundary condition at the aggregate surface of the sphere. In [|37| two other 
approaches were considered, a phenomenological model for the boundary 
condition and a radiation boundary condition. Both approaches allowed for 
incorporation of particle detachment in Smoluchowski model. Explicit ex- 
pressions for the concentration and intake rate of particles were given in the 



long time limit [37 



In this paper we consider a single, motionless three-dimensional droplet 
growing by adsorption of diffusing monomers on a two-dimensional substrate. 
The diffusing non-coalescing monomers are adsorbed at the aggregate perime- 
ter of the droplet with different boundary conditions. Models with different 
boundary conditions for the concentration of monomers are considered and 
solved in a quasistatic approximation. For each model, the diffusion equa- 
tion is solved exactly, subject to a fixed boundary. Using mass conservation 
law at the aggregate perimeter of the growing droplet, we then obtain an 
expression for the growth rate of the moving boundary. Explicit asymptotic 
solutions in the both short and long time limits are given for the concen- 
tration, total flux of monomers at the perimeter of the growing droplet and 
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for the growth rate of the droplet radius. This paper is organised as follows. 
In section 2, a model with an adsorption boundary condition is examined. 
In sections 3 and 4 we consider the two approaches which were introduced 
in [[37j to allow for particle detachment. A phenomenological model and a 
model with a radiation boundary condition are considered in sections 3 and 
4, respectively. Another boundary condition which assumes a constant flux 
of monomers at the aggregate perimeter of the droplet, is also introduced in 
section 5. Finally, in section 6 we compare the results of different approaches 
and summarise our conclusions. 

2 Growth Equations with Adsorption Bound- 
ary Condition 

Consider an immobile three-dimensional droplet which is initially surrounded 
by monodisperse droplets. The droplet lies on a two-dimensional plane sub- 
strate on which the monomers diffuse. Monomers have the volume V and 
diffuse with the diffusion constant D. Then, the concentration of monomers 
at point r and at time t, c(r,t), is described by the diffusion equation 



for r > R, where R(t) is the radius of the immobile growing droplet. The 
initial conditions are given by 



dc(r, t) 
dt 




(1) 




(2) 



which is the initial, uniform, monomer concentration and 
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R(t = 0) = 0, (3) 

which shows that the droplet is not present at the beginning of the process. 
We consider an adsorption boundary condition at the perimeter of the droplet 



c(r = R,t > 0) = (4) 

and assume that at infinity the concentration of the monomers is finite and 
equal to Cq. Concentration gradients in the neighborhood of the droplet 
create a flux of monomers on the two-dimensional substrate. This flux feeds 
the growth of the droplet. Therefore, the rate of increase of the droplet 
volume is related to the total flux of monomers at the perimeter of the droplet 
by mass conservation, 



$(t) = \R 



dR 



(5) 



where the total flux 



= V 



dc 

2tiRD — 
or 



R. 



(6) 



corresponds to the monomers incorporated at the perimeter of the droplet. 
In (|5D A is a dimensionless factor related to the contact angle of the droplet. 

In order to solve ([I]) with (2-4), we introduce the Laplace transform of 
the concentration, 



dt e _st c(r,t), 



(7) 



which satisfies the equation 



1 d ( dc\ _ 
r or \ or I 



(8) 
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Here we have already used the initial condition (^). The general solution of 
this equation is given by 



c(r, s) 



c 



+ A{s)K {qr) + B(s)I (qr), 



(9) 



where q = y s/D, and K and Iq are Modified Bessel functions of order zero. 
To have a finite solution as r — > oo, we set -B(s) = 0. The boundary condition 
@) in the Laplace transform version becomes 



c(r = R,s) = 0. (10) 

Using ([ToD the transformed concentration and its gradient normal to the 
droplet perimeter, yield 



c(r, s) 



dc(r, s) 



Co 
s 



1 - 



Co 



Kp{qr) 
K (qR) 

Ki{qr) 



'121 



9r (Ds) 1 / 2 K {qR)' 
To find time dependent concentration and its radial gradient, we use the 
Inversion theorem for (|IT] , |T2"|) . Both (|IT] , |T2"|) have a branch point at s = 0, so 
in the Inversion formula, we use a contour which does not contain any zeros 
of s and K (qR). Consequently, time dependent concentration and also total 
flux at the droplet perimeter from (Rf), are given by 



2co 

7T JO 



r(r.l) ^ I * < 1 



J {Ru)N (ru) - J (ru)JV (i2u)l rfw 



(13) 



*(*) 



8coW 



7T 



du 



[J$(Ru) + Ng(Ru)] u 



(14) 
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where J and N are Bessel functions of order zero. Using (|§,0) a differential 
equation for the growth rate of the droplet radius can be obtained 



XR 



dR 8c DV 



-Du 2 t 



du 



(15) 



dt 7T Jo [Jq(Ru) + Nq(Ru)] u ' 

which gives a general solution for R as a function of the time. We are 
interested in the short and long time solutions for the concentration, the 
total flux of monomers at the perimeter of the droplet and the growth rate 
of the droplet radius. 

For small values of the time, it is shown that the behaviours of c(r, t) and 
dc(r,t)/dr may be determined from the behaviors of c(r,s) and dc(r,s)/dr, 
respectively, for large values of the transformed parameter s. Then, we ex- 



pand the Bessel functions occuring in (fTLyi2 ) supposing s to be large. The 
final result for the concentration of monomers, keeping the leading time de- 
pendent term, is 



c(r, t) ~ c 



rj \2sjDt) 



(16) 



The total flux at the droplet perimeter and the growth rate of the droplet 
radius also in this limit using (|5|) and (pp, respectively, are given by 



$(t) ~ 2 c VRV^ t~ 1/2 , (17) 




We see that in the short time limit, R grows as a power of the time with an 
exponent of 1/4. 

For large values of the time, the behaviours of c(r, t) and dc(r, t) / dr may 
be determined from the bahaviours of c(r,t) and dc(r,s)/dr, respectively, 
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for small values of the transform parameter s. We then expand the Bessel 
functions occuring in (|TT|JI^) supposing s to be small. Keeping the leading 
time dependence term, the concentration of monomers yields 



In 



c(r, t) ~ 2 c 



R 



In 



ADt 
'a 2 ~R? 



(19) 



where a 



1.78107..., where 7 = 0.57722... is Euler's constant. The 



total flux at the droplet perimeter and the growth rate of the droplet radius 
also in this limit using (g) and (0), respectively are given by 



$(t) ~ AucqDV 



R(t) ~ A 



In 



4L>i 



1/3 



-1 



ln(r) 



(20) 



(21) 



where A = (9irVa 2 / A) 1/3 and r = 4c Dt/a 2 is the dimensionless time. Up 
to a constant, these are the same results which were obtained by Krapvisky 
based on a quasistatic approximation using a similarity variable approach 

PL 



3 Phenomenological Rate Equation Model 

One can consider various modification of the initial and boundary conditions 
(2) and (4). Here we improve the model and incorporate effects other than the 
irreversible adsorption at r = R expressed by (4), Ref. [07[| . In this section 
we consider a phenomenological modification of the boundary condition (4) 
to allow for detachment. This was introduced in |T7[] where the relation 
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?E^l = - mc (r } t) + k (22) 

at r = R was replaced for (4) . Here it is assumed that the diffusing monomers 
that reach the perimeter of the droplet are incorporated in the aggregate 
structure at the rate mc proportional to their concentration at R. The second 
term in (|22] ) corresponds to detachment and is assumed that only depends 
on the internal processes, so there is no dependence on the external diffuser 
concentration |3j 

To solve (1) with (2,3) and (0), we go through steps similar to section 2 
and only emphasize the final expressions. In the Laplace transform version, 
the boundary condition becomes 

k 

(s + m) c(r, s) = - + c (23) 
s 

at r = R. The concentration and the radial gradient of the concentration in 
this version become 

c mc -k K (qr) 

c(r,s) = (24) 

s s(s + m) K {qR) 

dc(r, s) _ mc - k 1 Ki(qr) 

dr ~ (Ds) 1 ' 2 (s + m) K (qR)' 1 ' 

Now we look for the solutions in the short and long time limits. 

For small values of the time, we use the asymptotic expansions of the 



Bessel functions in ( P4| , |25[) for large values of s and ignore m in comparison 
to s in the term (s + m). Then, the concentration, the total flux at the 
droplet perimeter and the growth rate of the droplet radius in this limit, 
keeping only the leading time- dependent terms, are given by 
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c(r, t) ~ c + 4mt f c - - £r/c ^= , (26) 



~ 4myi2V7rD c It 



1/2 



(27) 



J2(f) 



16mV , — — / k 
—yfrDU-- 



1/2 



(28) 



We see that in a phenomenological model, R grows as a power of the time 
with an exponent of 3/4 in the short time limit. In the expressions (|26| - p8D , in 
comparison with ([16| - |T8D in the previous section, there is a term as (co — k/m) 
which shows a reduction of the rate due to detachment, proportional to the 
ratio k/m. 

For large values of the time, we use the expansions of the Bessel functions 
in (^4|, ^5|) supposing s to be small and ignore s in comparison with m in the 
term (s+m). Then, the concentration, total flux at the droplet perimeter and 
growth rate of the droplet radius, keeping only the leading time dependent 
terms, yield 



c(r,t)~ - + 2(c -- 
m \ m , 



In 



$(t) ~ AnDV I c - — I fin 
V m j 



ADt 
~& 2 R 2 . 

ADt " 1 



a 2 R 2 



(29) 



(30) 



R(t) ~ A 



11/3 



ln(r) 



(31) 



where A = (9tcVct 2 / X) 1 / 3 and r = ADt(c — k/m) /a 2 . These asymptotic 
expressions are quite similar to the long time forms (|i~9l - pi~l) in section 2. The 
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only change is the reduction of the rate due to the detachment, proportional 
to the ratio k/m. 

For a fast enough detachment, the total fluxes of the monomers at the 
boundary in the both short and long time limits (|27| ,30) can actually become 
negative. In this case, the flux does not feed the growth of the droplet and the 
droplet volume does not increase anymore. Therefore, the mass conservation 



@ does not hold and the growth laws (|28| , |3T|) are not valid anymore. For a 
case in which cq = k/m, the system reaches a stationary state and therefore 
the total rate and the total flux of the monomers at the droplet perimeter, 
become zero for all times. Consequently, there is no growth for the droplet 
and the concentration of the monomers is equal to the initial concentration, 
Co, for all times. These results can be obtained from the both short and long 
time expressions (MM) and <MM), respectively. 



4 Radiation Boundary Condition 

In this section we consider another modification of the boundary condition 
(4) and replace it with a radiation boundary condition 

a ^M +/3 = c(M) (32) 



at r — R, Ref . |37[ . Here it is assumed that the concentration is proportional 
to its derivative, with an additional constant (3. Again we go through steps 
similar to the section 2 and only emphasize the final expressions. In the 
Laplace transform version, the boundary condition becomes 

dc{r, s) 8 

a ~fr~ + - = c (^ s ) ( 33 ) 
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at r = R. The concentration and its radial gradient in this version become 



Co _ (cp - g) -Kp^r) /^v 

l ' J s s K (qR) + aqKx(qR) ' 1 J 

ggr, a) = (cp - g) -K~i(gr) 

cV {Ds) 1 ! 2 K {qR) + aqK^qR)' 1 J 

We concentrate our attention to the solutions in the short and long time 

limits. 

For small values of the time, we use the asymptotic expansions of the 
Bessel functions in to get the leading time dependent terms for the 

concentration, the total flux and the droplet growth rate 

/x 2i, n ,/DRt\ 1/2 ^ „ fr-R\ 
c( r ,t)~c --{c -/3)( Erfcl-^), (36) 



a \ r J \2y/Dt) 



. ^ 2ttRDV . > , , 

m ~ (c -/?), (37) 



4ttW 



(co - /?) 



1/2 



4/2 



(38) 



aX 

The term (co — (3) in these expressions shows a reduction of the rate due 
to the detachment, proportional to the ratio (3. We see that in the short 
time limit, the total flux at the droplet perimeter is time-independent and R 
grows as a power law with an exponent equal to 1/2. 

For large values of the time, we expand the Bessel functions in ( |5^J33D 



supposing s to be small. Consequently, the asymptotic expressions for the 
concentration, total flux and the droplet growth rate, yield 



c(r, t)~/3 + 2(c - 13) - \% - , (39) 
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In 


KRJ 




In 




Ka 2 R 2 J 



$(t) ~ AttDV(co - p) 



In 



4L>t 



(40) 



ln(r) 



1/3 



(41) 



where A = (9irVa 2 / X) 1 ^ 3 and r = 4Dt(co — (3)/<J 2 - These long time expres- 
sions have the same forms as (p^ - [3lD provided we identify 



k 
m 



(42) 



For a fast enough detachment, analogue to the section 3, the total fluxes of 
the monomers at the boundary ( p7|J40| ) can become negative. In this case, 
the growth laws (|3"8|j4~T|) do not hold anymore. For a case in which cq = (3, 
analogue to the section 3, the system reaches a stationary state and therefore 
the total flux and the droplet growth rate become zero. The concentration 
also does not change with the time and is equal to the initial one. These 
can be seen from the both short and long time results (|36| - |3"8D and (|39| - ^TD , 
respectively. 



5 Constant Flux Boundary Condition 

In this section we impose a condition on the flux of the monomers assum- 
ing that the total flux of monomers at the droplet perimeter is constant. 
Therefore, we replace (4) with 

$(t) = Q (43) 



13 



at r = R, where $(t) is given by ([]) and Q is a constant. The analogue to the 
previous sections, in the Laplace transform version, the boundary condition 
becomes 



2nRDV^^ = ® (44) 
or s 

at r = R. The concentration and its radial gradient in this version are 

c( r s) = ^ 9. (45) 

and 

dc(r, s) = Q Kxjqr) 

dr 2-kRDV sK 1 (qR)' 1 ' 

Appropriate expansions of the Bessel functions in (^5j) give us the limiting 

forms of the concentrations. For small values of the time it yields 

c ( r , t ) ~ c - % -% ( —) Erfc ( (47) 
and for large values of the time it gives 

-JwH^)- (48) 

The trivial solution for the droplet growth rate using (|],f|3]) is 

R(t)=( 3 4-) 1/3 t 1/3 i (49) 



A J 



for all times. 
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6 Conclusions 



We studied the growth of a single, motionless, three-dimensional droplet that 
accommodates monomers at its perimeter on a 2D substrate. The noncoa- 
lescing monomers diffuse and are adsorbed at the aggregate perimeter of the 
droplet with different boundary conditions. Models with adsorption and radi- 
ation boundary conditions, and a phenomenological model for the boundary 
condition, were considered and solved in a quasistatic approximation. In a 
model with adsorption boundary condition, the droplet forms an absorber 
and the concentration of the monomers at its perimeter is zero. In a phe- 
nomenological model, we assumed that the diffusing monomers that reach 
the perimeter of the droplet, are incorporated in the aggregate structure at 
a rate proportional to their concentration at the boundary. We also added 
another term which corresponds to detachment. In a model with radiation 
boundary condition we assumed that the concentration is proportional to its 
derivative with an extra detachment term. For each model, we solved ex- 
actly the diffusion equation for the concentration of the monomers, subject 
to a fixed boundary. Then, using a mass conservation law at the perimeter 
of the droplet, we found an expression for the growth rate of the moving 
boundary. Models were subjected to an initial, uniform concentration of 
monomers. Asymptotic results for the concentration, total flux of monomers 
at the boundary and the growth rate of the droplet radius, were obtained 
in both short and long time limits. The results revealed that in both phe- 
nomenological and radiation models, in comparison with adsorption model, 
there is a reduction of the rate due to the detachment. The rate can become 
negative if the detachment is fast enough. In this case, the total flux of the 
monomers at the perimeter of the droplet become negative. Therefore, the 
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flux does not feed the growth of the droplet volume and the droplet growth 
laws obtained in the models, are not valid anymore. For a value of the de- 
tachment for which the total rate and therefore the total flux become zero, 
the system reaches a stationary state and there is no growth for the droplet 
anymore. The same reduction of the rate was obtained in [37| where incor- 
poration of particle detachment in Smoluchowski model of colloidal growth, 
was considered. 

The results in the short time limit predicted that the radius of the droplet 
grows as a power of the time with different exponents for different bound- 
ary conditions. The exponents of the power laws were 1/4, 1/2 and 3/4, 
respectively, for the models with adsorption, radiation and phenomenologi- 
cal boundary conditions. We see that the growth rate is the slowest for the 
adsorption boundary condition and is the fastest for the phenomenological 
model. This is because, as was said before, in the phenomenological model, 
the diffusing monomers at the perimeter of the droplet are incorporated in 
the aggregate structure of the droplet. The total flux of the monomers at the 
droplet perimeter is also power law with exponents of —1/2 and 1/2 for the 
adsorption and phenomenological model, respectively, and is a constant for 
the radiation model. Again the flux is maximum for the phenomenological 
model and is minimum for the adsorption model. 

In the long time limit, the growth law for the radius of the droplet was 
the same for all boundary conditions. Also the concentration and total flux 
had the same time dependency in all models. The only change, as we said 
before, was the reduction of the rate due to the detachment in the both phe- 
nomenological and radiation models in comparison with adsorption model. 
Asymptotic results for large values of the time exhibited that the radius 
of the droplet increases as [t/ ln(t)] 1 / 3 in all models. This was obtained by 
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Krapivsky |[TjJ where a similarity variable approach was used to treat the 



growth of a droplet with an adsorption boundary condition based on a qua- 
sistatic approximation. 

We saw that the time dependency of the results was the same for all the 
models in the long time limit and was different for different models in the 
short time limit. This suggests that initially the flux of the monomers at the 
boundary and therefore the droplet growth rate, are affected by the condi- 
tion at the boundary. But in the long time limit, the system reaches a stable 
state and the initial effects can be ignored, therefore all the models give the 
same results. This suggests that a rate as [t/ln^)] 1 / 3 is a universal asymp- 
totic growth law for the radius of the droplet independent of the boundary 
conditions. 

In the both models with phenomenological and radiation boundary con- 
ditions, similar to the results in [[37]], the value of the concentration at r = R 
for large times, see (^9|) and (|39"D, is exactly equal to (3 = k/m, independent 
of R. This suggests that as far as large R and large time behaviours are 
concerned, we can use the boundary condition 

c = /3=- (50) 
m 

at r — R, instead of phenomenological and radiation boundary conditions. 
Indeed, the value of the concentration at r = R is the only parameter needed 
to calculate the modifications of the asymptotic behaviours due to the de- 
tachment. With this boundary condition for all times, the asymptotic results 
(]29| - |3~T|) and ( p9| - |4l|) become exact. Therefore, a constant concentration of the 



monomers at r = R for all times, gives an exact growth rate as [t/ ln(t)] 1 / 3 . 

We also examined another model with a constant flux of monomers at 
r = R. The results showed that the radius of the droplet grows as t 1 / 3 for all 
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times. Thus, the growth laws predicted by a constant concentration and by 
a constant flux at the boundary, differ from each other by a slowly varying 
logarithmic factor. 
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